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ABSTRACT 

We consider the luminous blue variable (LBV) star V532 in M33, also known as 
Romano's star, in two different spectral states: in the optical minimum of 2007/2008 
and during a local brightening in 2005. Optical spectra of low and moderate resolution 
are modeled using the non-LTE model atmosphere code CMFGEN. All the observed 
properties of the object in the minimum are well described by a late nitrogen-sequence 
Wolf-Rayet (WN) star model with relatively high hydrogen abundance (H/He= 1.9), 
while the spectrum during the outburst corresponds to the spectral class WNll and 
is similar to the spectrum of P Cyg. The atmosphere is enriched in nitrogen by about 
a factor of 6 in both states. Most of the heavy element abundances are consistent with 
the chemical composition of M33. Bolometric luminosity is shown to vary between the 
two states by a factor of ^1.5. This makes V532 another example of an LBV that 
shows variations in its bolometric luminosity during an outburst. 

Key words: galaxies: individual: M33 - stars: Wolf-Rayet - stars: supergiants - 
stars: individual: Romano's star (M33) 



1 INTRODUCTION 

Luminous blue variables (LBVs) are rare objects of very 
high luminosity (~ l(f Lq) and mass loss rates (typically 
10~®MQyr~^, up to 10~'*MQyr~^), exhibiting strong ir- 
regular photometric and spectral variability (Conti 1984; 
Humphreys & Davidson 1994). They are generally believed 
to be a relatively short evolutionary stage in the life of a mas- 
sive star, marking the transition from the Main Sequence to- 
ward Wolf-Rayet (WR) stars (Humphreys & Davidson 1994; 
Smith, Crowther & Prinja 1994; Maeder & Meynet 2000). 
Significant part of the stellar mass is lost during this stage 
forming massive circumstellar nebulae (Nota et al. 1995; 
Weis 2001) like the spectacular Homunculus surrounding 
r] Car (Smith et al. 2003). However, recent investigations 
of the light curves of a few supernovae (SN2006gy, 2006tf, 
SN2005gl) indicate that their progenitors underwent LBV- 
like eruptions (Smith & McCray 2007; Smith et al. 2008; 
Gal- Yam & Leonard 2009). These observations support the 
view that at least some luminous LBV stars are the end 
point of the evolution and not a transition phase, but which 



* Based in part on data collected at Subaru Telescope and ob- 
tained from the SMOKA, which is operated by the Astronomy 
Data Center, National Astronomical Observatory of Japan, while 
other data were taken from the archive of Special Astrophysical 
Observatory (SAO) of Russian Academy of Sciences (RAS) 
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is still in contradiction to current evolution models (Meynet 
et al. 2011). 

LBVs span quite a large range of magnitudes and vari- 
ability types (Humphreys & Davidson 1994). Some of them 
show "normal eruptions" also called an S Dor variability 
phase. During this phase the star may have brightness of 
up to 2 mag, but the bolometric luminosity remains essen- 
tially constant. Other LBVs occasionally experience "giant 
eruptions" , during which they increase their luminosity and 
can reach M « — 14mag (4-6 mag brighter than their typi- 
cal quiescence magnitudes) . Of the 35 Galactic LBVs (con- 
firmed and candidate) (Clark et al. 2005) only three {rj Car, 
P Cyg and AFGL2298) were observed during giant eruptions 
(Humphreys & Davidson 1994; Clark et al. 2009). Several ex- 
amples of eruptions accompanied with changes in bolomet- 
ric luminosity have been studied spectrally (such as HD5980 
(Koenigsberger 2004), NGC 2363-Vl (Drissen et al. 2001), 
AFGL2298 (Clark et al. 2009)). These giant eruptions are 
not only spectacular events but are probably responsible for 
bulk of the mass loss by very massive stars (Smith & Owocki 
2006). 

Observations of Galactic LBV stars are inevitably con- 
nected with difficulties in determination of the distance and 
interstellar extinction, that results in huge uncertainties in 
their bolometric luminosities. Hence, studying these rare 
objects in nearby galaxies is potentially more prospective, 
though extragalactic objects are more distant and hence 
more difficult targets for spectral observations. 

Romano's star is located in the outer spiral arm of 
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Figure 1. The B-band light curve of V532. Arrows indicate the 
dates of spectral observations. 



the M33 galaxy at a distance of about 17' from the cen- 
tre. The first light curve for V532 was presented by Romano 
(1978), while the first spectral observations were obtained 
in 1992 (Szeifert 1996). Because V532 demonstrates pro- 
nounced photometrical and spectral variability (Kurtev et 
al. 2001; Polcaro et al. 2003), it is classified as an LBV. 
The object changes from a B emission line supergiant in 
the optical maximum (Szeifert 1996), through Ofpe/WN 
(WN10,WN11) and WN9 toward a WN8-like spectrum in 
deep minima (Maryeva & Abolmasov 2010). Figure 1 shows 
a B-band light curve of V532 from 1990 to 2008. 

During 2004 and 2005, V532 becomes brighter by about 
1™ in B and reaches 17™ in this band (see Zharova et al. 
(2010), Maryeva & Abolmasov (2010), and Polcaro et al. 
(2010) for details). We classify the spectra obtained dur- 
ing this period as WNll, in agreement with the estimates 
by Polcaro et al. (2010). Starting from the middle 2005, 
Romano's star weakens in all bands. Its visible magnitude 
reaches 18™. 68 in the V band in February 2008. Recently, it 
exhibits a slight brightening by about 0™. 2 (Polcaro et al. 
2010). We estimated the spectral class of early 2008 to be 
WN8 (Maryeva & Abolmasov 2010). 

In this work we investigate the optical spectra of V532 
in two different states, the brightness minimum of 2008 and 
a moderate brightening in 2005, using the non-local ther- 
modynamic equilibrium (non-LTE) radiative transfer code 
CMFGEN (Hillier & Miller 1998). We refer to both states as 
hot and cool phases, respectively. 

This paper is organized as follows. Observational data 
and data reduction process are described in Section 2. In Sec- 
tion 3 we describe the basic properties of the cmfgen code. 
We devote Section 3.2 to characteristic diagrams, while in 
Section 3.3 and Section 3.4 we present and analyse the mod- 
eling results for hot and cool phases of V532. In Section 4 
we discuss the results. Finally, Section 5 we summarize the 
main points of our work. 



2 OBSERVATIONS AND DATA REDUCTION 

In this work, we use one spectrum of Romano's star ob- 
tained during the outburst in 2005 at the Special Astro- 
physical Observatory (SAO) 6-m telescope ^. For studying 
Romano's star at minimum brightness we use one spectrum 
obtained in .lanuary 2008 at the 6-m SAO telescope and 
one spectrum obtained at the SUBARU telescope ^. The 6- 
m telescope data were obtained with the Spectral Camera 
with Optical Reducer for Photometric and Interferometric 
(SCORPIO) in the long-slit mode (Afanasiev & Moiseev 
2005). One exposure 1200 s in length was obtained with 
the SUBARU telescope with the Faint Object Camera and 
Spectrograph (FOCAS) (Kashikawa et al. 2002) in October 
2007. VPHG450 grism was used providing the spectral range 
of 3750-5250 A. Slit width of 0!'5 implies spectral resolution 
of about 1.7 A . Observational log information on the data 
used in this work is summarized in Table 1. 

All the SCORPIO spectra were reduced using the ScoRe 
package. ScoRe was written by Maryeva and Abolmasov 
in IDL language for SCORPIO long-slit data reduction. 
This package consists of procedures written by V. Afanasiev, 
A. Moiseev, P.Abolmasov and O. Maryeva. Package includes 
all the standard stages of long-slit data reduction process. 
FOCAS data were reduced in IDL development environment 
using procedures similar to those consisting ScoRe but tak- 
ing into account the specific features of FOCAS. We describe 
the observational data and data reduction process in more 
detail in Maryeva & Abolmasov (2010). 

3 THE MODEL 
3.1 CMFGEN code 

For our analysis we used the non-LTE radiative transfer 
code CMFGEN (Hillier & Miller 1998). CMFGEN has been ap- 
plied to several classes of objects where non-LTE effects and 
stellar wind are important (e.g., WR, LBV, and O stars). 
More recently, CMFGEN was used to investigate the pho- 
tospheric phase of Type II supernovae (Dessart & Hillier 
2005). CMFGEN solves radiative transfer equation for objects 
with spherically-symmetric extended outflows using either 
the Sobolev apijroximation or the full comoving-frame solu- 
tion of the radiative transfer equation. To facilitate simulta- 
neous solution of the transfer equation and statistical equi- 
librium equations, partial linearization method is used. To 
facilitate the inclusion of metal line blanketing in CMFGEN, 
supcrlevcl approach (Anderson 1985, 1991) is used. In this 
formalism, levels with similar properties are treated as one 
and have identical departure coefficients. This allows to save 
a considerable amount of computer memory and time. Re- 
cent versions of the code incorporate also the effect of level 
dissolution, influence of resonances on the photoionization 
cross section, and the effect of Auger ionization. 

Clumping is incorporated into CMFGEN using volume 

^ Spectral data were taken from the aicliive of Special Astrophys- 
ical Observatory (SAO) of Russian Academy of Sciences (RAS) 
(http : //www. sao .ru) 

^ Spectral data from the SUBARU telescope were taken from the 
SMOKA science archive (Baba et al. 2002) (http://smoka.nao. 
ac.jp) 
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Table 1. Observational log for the spectral data used in the work. S/N is signal-to-noise ratio per resolution element. 
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Exposure Spectral Spectral 

Date time Grism range S\ S/N Seeing standard 

[s] [A] [A] ["] star 



06.02.2005 SCORPIO 2 x 300 VPHG550G 3700-7200 10 8 1.7 G248 

08.01.2008 SCORPIO 2 x 900 VPHG1200R 5700-7500 5 20 2.1 BD25d4655 
08.10.2007 FOCAS 1200 VPHG450 3750-5250 1.7 30 0.53 BD40d4032 



filling factor approach (Hillier & Miller 1999). Filling fac- 
tor is allowed to depend on radius. By default, the wind is 
considered homogeneous at the hydrostatic radius and bo- 
comes more and more clumped with the wind velocity. Tak- 
ing clumping into account decreases the derived mass loss 
rates by a factor of ~ 3 — 5. The unclumpod and clumped 
mass loss rates are related to the volume filling factor / as 
Muncl = Mci/^/f. 

Each model is defined by the hydrostatic stellar radius 
J?t, luminosity L*, mass-loss rate Mci, filling factor /, wind 
terminal velocity Voo, stellar mass M, and by the abun- 
dances Zi of included species. Hydrodynamic equations are 
not solved, instead we propose to use /3-law for velocity and 
fix mass-loss rate for each model. This allows to define the 
density structure throughout the wind. 

3.2 Diagnostic Diagrams 

In order to reproduce the spectrum of V532 obtained in Oc- 
tober 2007 while the object was in a deep minimum we took 
the model of the WN8 star WR40 (HD96548) calculated 
by Herald ot al. (2001) and gradually changed the parame- 
ters of the model. We calculated a grid of about 130 models 
with difTerent parameter values (luminosity, mass-loss rate, 
mass, elementary abundances). Luminosity was varied in the 
ranee (0.4 ^ 2) X W'^Lq. For every model, we recalculated 
the model flux for the distance of M33. For our calcula- 
tions we adopted a distance to M33 of -D = 847 ± 60 kpc, 
which gives a distance modulus of (m — M) = 24.64 ± 0.15™ 
(Galleti et al. 2004). Then, the simulated spectra were con- 
volved with B- and V-band sensitivity filters. The resulting 
fiuxes were converted to magnitudes (Leng 1974) and com- 
pared to the photometrical data {B = 18™. 5 ± 0.05 and 
V = 18™. 68 ± 0.05). According to the detailed dust maps 
of M33 by Hippelein et al. (2003) the effect of the intrinsic 
extinction should be negligible near the object. Extinction 
in the Galaxy is estimated as E{B - V) = 0™. 052 by the 
NED extinction calculator (Schlegel et al. 1998). This value 
is well inside the errors for the observed B — V colour index 
but affects the observed luminosity values that should be 
thus increased by ~ 0.17 mag. 

Mass loss rate was varied in the limits (1.2 — > 7.7) x 
lO-'^Moyr-i typical for WN8-WN11 stars (Crowther et al. 
1995c; Hamann et al. 2006). In all the models we assume that 
the volume filling factor at infinity equals foo = 0.1. Since 
V532 resides in low-metallicity environment, we considered 
only sub-solar iron abundances (Fe/Fe© = 0.4 ± 0.15). We 
also expected sub-solar abundances for Si, Mg, Na. 



The velocity law used was assumed a simple /3-law with 
/3 = 1. Photospheric velocity was set to 100 km/s and the 
terminal velocity is 400 km/s for all the models used for 
diagnostic diagrams. Below we investigate the effect of vari- 
ations in /3 parameter and photospheric velocity on model 
spectra. Profile fitting for the triplet lines of He I (such as 
A3889, 4025, 4471) allows to estimate the terminal velocity 
as ~ 400 km s"^ (Maryeva & Abolmasov 2010). 

Each model was classified using the equivalent width 
(EW) ratio of emission components of HeiA5876 and 
HeiiA5411 (Smith, Shara & Moffat 1996). More precise es- 
timates of physical conditions in the wind and atmosphere 
may be made using EW ratios which weakly depend on spec- 
tral resolution and elementary abundances and also allow to 
exclude the contribution of a possible circumstellar nebula. 

We construct several characteristic diagrams to com- 
pare the model spectra with the observations using equiva- 
lent width ratios. Figure 2 shows one such diagram. Model- 
ing is complicated by the nebula surrounding V532. There- 
fore we use the characteristic EW ratios of He iiA4686 to 
HelA5876 (these lines form both in the stellar atmosphere 
and in the nebula) and HeiiA5411/HeiA4713 (contribution 
of the nebula to these two lines should be negligibly small). 
In this diagram, models fall along a narrow curve with the 
spectral class changing smoothly along it. Effective temper- 
ature changes from T* < 25kK for the lower left corner to- 
ward T, > 40kK for the upper right. Model distribution in 
a narrow locus along a single curve is probably connected to 
the correlated behavior of emission lines of He II as well as 
He I. Most models shown in the diagram have identical wind 
velocity and structure but different radii, luminosities and 
mass loss rates. We see that some real stars lie near this lo- 
cus while WR 108 is offset. WR 108 is an unusual WN9 star 
it has a higher terminal velocity (~ 1200km/s Hamann 
et al. (2006)) than other WN9 stars. Crowther et al. (1995a) 
classified WR 108 as intermediate star between normal Of 
and WN stars. 

Bright hydrogen lines are present in the spectra, ev- 
idently stronger than in the spectra of ordinary late WN 
stars. In this sense, Romano's star is similar to H-ricli WN 
stars. Therefore we calculate models with H/He= 0.75 — > 2.6 
(by number) that is typical for hydrogen WR stars. For the 
models that reasonably fit the observed equivalent widths, 
we varied the chemical composition and filling factor to fit 
the line strengths of individual elements and line profiles, 
correspondingly. 
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Figure 2. The plot of the EW ratio logarithm log(HeII 4686 /Hel 5876) versus log(HeII 5411 /Hel 4713). V532 location in October 
2007 (hot phase) is marked by a diamond. Sk-66 40 (WNIO), BE381 (WN9), HDE269927c (WN9) and WR108 (WN9-abs) are shown for 
comparison. Data on these objects were taken from Crowther et al. (1995b). 



3.3 Hot-Phase Spectrum 

In Figure 3 we show the observed spectra of V532 at mini- 
mum brightness and the best-fit model spectrum. Figure 3 
shows good agreement between the observed and model 
spectra. Both in observational data and in the best-fit 
model, singlet and triplet lines of Hex show different pro- 
file shapes: triplet lines have typical P Cyg profiles while 
singlets are pure flat-topped emissions. 

For this model, the best-fit volume filling factor at in- 
finity foo = 0.1, mass loss rate Mcl = 1-9 x lO'^^Moyr"^ 
(which correspond to unclumped Munci — 6x lO~^M0yr~^). 
The stellar temperature T* follows from the relation L, = 
AnR^aTf, and the effective temperature at the photosphere 
Teff is defined by the Rosseland optical depth of 2/3. The 
values obtained for T,, Te// and other physical parameters 
are listed in Tables 2 and 4. Errors in luminosity given in the 
table are determined by the uncertainties of the photomet- 
rical data. The errors in the H/He ratio are from the model 
fitting alone. 

Surface chemical abundances obtained for V532 are 
listed in Table 3. Abundance pattern is consistent with the 
moderately sub-solar metallicity of M33 ([Fe/H]~ —0.5), 
but nitrogen is significantly over-abundant (~ 6.4^© ). 
The latter value is consistent with the existing evolution- 
ary models and with the data on other nitrogen-rich WR 
stars (Herald et al. 2001). Silicon abundance was adjusted 
using Sim A4565.05 and Siiv AA4088.90, 4116.10 lines, mag- 
nesium abundance using Mgii A4481.13. Ne, Al, S, Ar and 
Ca abundances were fixed relative to He because they are 
poorly constrained by observational data due to the lack of 
strong lines for these elements. We assume these abundances 
identical to those for WR40. 

After we obtained a model reasonably approximating 



Table 3. Derived stellar abundances for Romano's star. Abun- 
dances for all the elements are identical for the two phases within 
the uncertainties. 



SPECIES 


Number Fraction 


Mass Fraction 


Zi/Z0 


H 


1.9 


3.18 X 10-1 


0.45 


He 


1.0 


6.7 X 10-1 


2.4 


C 


1.0 X 10-4 


2 X 10-4 


0.07 


N 


3.0 X 10-3 


7 X 10-3 


6.4 





8.0 X 10-" 


2.2 X 10-3 


0.23 


Ne 


2.4 X 10-* 


8.1 X 10-4 


0.47 


Na 


7.0 X 10-*^ 


2.7 X 10-^ 


0.78 


Mg 


5.0 X 10-^ 


2 X 10-4 


0.31 


Al 


8.4 X 10-6 


3.8 X IQ-"' 


0.7 


Si 


2.0 X 10-4 


1.0 X 10-3 


1.43 


S 


2.3 X 10-^ 


1.2 X 10-4 


0.34 


Ar 


6.5 X 10-« 


4.4 X 10-5 


0.43 


Ca 


4.2 X 10-6 


2.8 X 10-5 


0.44 


Fe 


4.87 X 10-5 


4.6 X 10-4 


0.34 



the observed spectra and consistent with the photometri- 
cal data we investigated the effects of photospheric veloc- 
ity on the model spectrum and the derived model param- 
eters. Decrease in photospheric velocity implies increase in 
Rosseland optical depth that affects both Te/f and J?2/3- 
These changes practically do not affect hydrogen and neu- 
tral helium lines, while the EW of the HeiiA4686 emission 
increases by about 14 per cent when photospheric velocity 
docroasos by a factor of two. Table 2 gives the parameters of 
a model with Vphot = W km s~i that reproduces the spec- 
trum indistinguishable from the current best-fit model with 
Vphot = 100 km s~i. One order of magnitude uncertainty in 
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Figure 3. Normalized hot-phase optical spectra compared with the best-fit CMFGEN model (dotted line). Bottom panel shows the 
spectrum obtained with SCORPIO, top and middle panels show the spectrum obtained with FOCAS. The model spectrum on bottom 
panel is convolved with the sA-wide Gaussian instrumental profile. 

Table 2. Best-fit model parameters for diff'erent values of /3 and photospheric velocity. H/He denotes hydrogen number fraction relative 
to helium. 



/3 


Vphot 




T* 


R, 




R2/3 






H/He 




[km/s] 


[km/s] 


[kK] 


[Ro] 


[kK] 


[R©] 


[105 Lq] 


[ lO""'M0yr-l] 




1 


100 


400 


34.0 ± 0.5 


20.8 ± 1. 


31.7 ±0.5 


23.9 ± 0.7 


5.2 ±0.2 


1.9 ±0.2 


1.9 ±0.2 


1 


10 


400 


33.0 ±1.0 


21.6 it 1.2 


29.4 ±0.5 


27.3 ±2.0 


5.0 ±0.25 


1.6 ±0.1 


1.9 ±0.2 


4 


100 


500 


33.0 ± 1.0 


21.6 ±1.2 


29.0 ± 1.0 


27.8 ±2.0 


5.0 ± 0.25 


1.6 ±0.15 


2.0 ±0.2 
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photospheric velocity results in 3 per cent and 10 per cent 
uncertainties in luminosity and mass loss rate, respectively. 



3.3.1 Hot-Phase Spectrum with /3 > 1 

As stated above the velocity law used was assumed a simple 

/3-law: 



with 



ro = R,{l 



/ Vphot ^ 



where Vphot is photospheric velocity. Above we fixed the /3 
parameter to unity. In some objects in certain phases, how- 
ever, (5 attains a significantly higher value, up to /3 ~ 4 
(Najarro et al. 1997; Groh et al. 2009). In general, LBV 
stars have low /3 < 1 during optical minima, but the veloc- 
ity law changes during eruptions (Guo & Li 2007). Let us 
now consider the effect of variable /3 parameter on model 
spectrum. 

Figure 4 illustrates the dependence of the 4630-4713 A 
blend shape on /3. Increasing /? produces higher equivalent 
widths of all the detectable lines of the blend. Besides, the 
wings of He iiA4686 (see also section 4.2) become stronger. 
As a result, intensities of C ill AA4647.4, 4650.3, 4651.5 Unes 
change, while the line ratios of C ill AA4647.4, 4650.3, 4651.5 
remain somewhat constant. For large /3 (/? ^ 2) ClIlA4647.4 
is stronger than other carbon lines, and the shape of the 
4630-4713A blend changes. Intensity ratios of the carbon 
lines observed in the spectrum of V532 (see the right-hand 
panel of Figure 4) indicate /3 ~ 4. An increase in /3 affects 
the model spectrum similarly as an increase in the mass- 
loss rate. Therefore we may describe the observed spectrum 
by a model with /3 = 4 and a lower (by about 10 percent) 
mass- loss rate. 

Figure 5 shows comparison between the observed spec- 
trum of Romano's star and the model with /3 = 4. The 
model with /3 = 4 fits better the triplet lines of helium, 
while He iiA4686 is brighter than for the best-fit model with 
P — 1. Table 2 shows the best-fit parameters for the models 
with different values of photospheric velocity Vphot and /3. 



3.4 Cool-Phase Spectrum 

During 2004 and 2005, V532 becomes brighter by about 

1 mag in B and reaches 17 mag in this band. Colour 
index B — V is constant within the observational errors 
{B - V = -0.11 mag for 2007 October and -0^.17 for 
February 2005). Spectra obtained during this brightening 
we classify as WNll, in agreement with the estimates of 
Polcaro ct al. (2010). For studying Romano's star in the 
cool phase, we use the spectrum obtained in 2005 February 
at the 6-m SAO telescope. 

Using the hot model as initial, we increased the mass 
loss rate and hydrostatic radius thereby decreasing the ef- 
fective temperature. After obtaining a model consistent 
with the photometrical data ( V = 17.27 ± 0.03mag, B = 
17.1 ± 0.03mag, assuming Ay = 0.17mag) we started vary- 
ing the wind velocity and volume filling factor. The EW 



of Ha increases by 30 per cent when fao is increased 
from 0.1 to 0.5. But the effective temperature at the pho- 
tosphere and the photosphere radius varies insignificantly. 
Mass loss increase leads to Ha EW increase as well. A 
18 per cent variation of mass loss rate leads to approxi- 
mately 18 per cent variation in the EW value. However, 
it also leads to significant changes in effective temperature 
and photosphere radius. Therefore we fixed M and changed 
only /oo- To compare the model spectrum with the obser- 
vations we convolved the observed spectrum with the lOA- 
wide Gaussian instrumental profile. Figures 6 and 7 show 
the observed spectrum and the best-fit model spectrum con- 
volved with the instrumental profile. The figures also show 
comparison between the raw (high-resolution) best-fit model 
spectrum and a spectrum of the LBV star P Cyg (Blla"*") 
obtained in August 1998 and taken from Elodie archive 
(http: //atlas . obs-hp.fr/elodie/ intro .html). The spec- 
tra arc unexpectedly similar. 

For the best-fit model in the cool phase, volume filling 
factor at infinity /oo =0.5. This value is factor of 5 higher 
than the one typically found for WR stars (Herald et al. 
2001). Note that this value is equal to that for P Cyg (Na- 
jarro 2001), thus confirming the similarity of these objects. 
The volume-filling factor seems to depend strongly either on 
the mass loss rate or on the spectral state, changing from 
~ 0.1 for late WN stars to ~ 0.5 for B-type hypergiants. 
Self-consistent modeling of wind acceleration is needed to 
understand the mechanisms leading to the strong dumpi- 
ness of the winds of Wolf-Rayet stars. 

The mass loss rate for this model is Md = (4.5 ± 
0.2) X lO-^Moyr-i, luminosity is (7.7 ± 0.25) x lO^L©, 
Te// = 20.4 ± 1.0 kK. Wind and stellar parameters of V532 
in the maximum of brightness are given in Table 4. Sur- 
face chemical abundances are the same as for the hot-phase 
model (see Table 3). 



4 DISCUSSION 

4.1 Evolution of the Physical Parameters 

We model the spectra of V532 in majcimum of bright- 
ness (V'=17™, Feb. 2005) and in minimum of brightness 
( ^=18™. 6, Oct. 2005) using the non-LTE radiative trans- 
fer code CMFGEN. Stellar parameters derived for both hot- 
and cool-phase models are given in Table 4. For compari- 
son, the values of these parameters for some other WN stars 
taken from the literature are given in the table. Note that 
the parameters of the Galactic WN8 stars WR124, WR16 
and WR40, as well as the LBV AG Car in photometric min- 
ima and LBV P Cyg were calculated using CMFGEN models 
taking clumping into account (see Crowther et al. (1999); 
Herald et al. (2001); Groh et al. (2009); Najarro (2001)). 
On the other hand modeling of the WN9h stars R84 and 
BE381 did not account for clumping, and that may lead 
to mass loss rate overestimates. Table 4 shows that in mini- 
mum brightness, V532 is similar to a classical WN8 star, but 
the wind velocity is lower, characteristic instead for a WN9 
star. We observe that relative hydrogen abundance (H/He) 
for V532 is similar to that of WN8h stars, as given in the 
table. V532 as well as other LBVs AG Car and P Cyg (Groh 
et al. (2009); Najarro (2001)) shows significant enhancement 
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Figure 4. Left-hand panel: effects of pajrameter /3 on the profile of the 4630-4713 A blend. Right-hand panel: fitting of the 4630-4713 A 
blend. Observational spectrum is shown by a solid line, and the dotted line shows our hot-phase model calculated with cmfgen (/3 = 4) . 



of helium (more than a factor of 2 relative to solar) , impying 
the end of the hydrogen shell-burning phase. High nitrogen 
content and depletion of carbon and oxygen are indicative 
of material which has undergone the CNO cycle. Our re- 
sults agree with the contemporary understanding of V532 
and Galactic LBV stars as a transitional stage from main- 
sequence supergiants to WR stars. 

During the 2005 February outburst, parameters of the 
star correspond to the spectral class WNll. The model spec- 
trum is similar to the spectrum of P Cyg. V532 shows a 
WNll spectrum in the maximum, while the classical LBVs 
such as AG Cax and P Cyg show the same spectrum only 
in the deep minima and the long-lasting quiet state, respec- 
tively. Note however that V532 had a strong maximum in 
1993 (0™. 9 brighter than in February 2005) and exhibited 
a B-supergiant spectrum. Hence, V532 shows stronger spec- 
tral variability than AG Car. 

The two phases, hot and cool, arc mainly distinguished 
by the photosphere radius. In the hot phase the radius is 
about three times larger in the cool phase. Three basic par 
rameters vary simultaneously and make measurable contri- 
butions to the observed inflation of the star. For the two 
states, AI differ by a factor 2.4, and the wind velocity is 1.8 
times larger for the hot state. It is easy to check that the 
size of the wind photosphere should scale approximately as: 



Rph oc 

' Voo 



where n = 1 in the case of pure scattering wind (mass- 
absorption coefficient x = const) and n = 2/3 when x oc p, 
as for true absorption processes. Hence, R^h is expected to 
vary by a factor 3 — > 5, in consistence with the observed 
value of 2.9. Our model favours correlation of hydrostatic 
radius with mass loss rate. 

Romano's star is situated near a young OB association 
OB 89. Probably V532 was a member of OB 89 and was 
ejected via slingshot-type dynamical interaction. More detail 
we consider this suggestion in Maryeva & Abolmasov (in 
preparation). 

Polcaro et al. (2010) presented the values of bolo- 



metric luminosity, effective temperature and radius of Ro- 
mano's star, using bolometric corrections for known WN8 
and WNll stars. Here, we use a more comprehensive way to 
estimate the bolometric luminosity. But the main conclusion 
holds: bolometric luminosities of V532 were different in 2005 
and 2008. The luminosity of V532 in 2005 (L. = 7.7- lO^'Lo) 
is 1.5 times higher. Therefore, V532 should be considered 
one more LBV (after the objects mentioned by Kocnigs- 
berger (2004); Drissen et al. (2001); Clark et al. (2009)) that 
changes its luminosity during (even moderate amplitude) 
eruption. In this sense, V532 behaves similarly to AG Car 
that has bolometric luminosity variations during its S Dor 
cycle (Groh et al. 2009). 



4.2 The Broad Wings of the HeiiA4686 Emission 

Polcaro et al. (2010) detect a broad (~ 1000km s~^) com- 
ponent of the HeIlA4686 line in the spectrum of V532 ob- 
tained in December 2008. Authors explain this component 
by a bimodal stellar wind. Although cmfgen adopts spher- 
ical symmetry, the model profiles of this emission line pos- 
sess similar, and sometimes even broader and brighter wings 
(Fig. 8). Observed wings may be explained by a more com- 
mon phenomenon - by electron scattering. Strong and very 
broad emission wings were first found for LBVs in P Cygni 
by Bcrnat & Lambert (1978). Balmer lines in the spectrum 
of AG Car do have such wings, which extend to more than 
±1500 km s"^ from the core of the line (Stahl et al. 2001). 
For P Cyg, electron scattering wings were reproduced by 
Najarro et al. (1997) in spherical symmetry. The wings are 
explained by scattering of line photons by free electrons of 
the stellar wind. Large widths of the wings reflect high ther- 
mal velocities of free electrons and does not correspond to 
bulk motion of the gas. 



5 CONCLUSIONS 

Using comoving frame numerical radiative transfer with the 
CMFGEN code, we estimate the physical parameters of the 
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Figure 5. The normalized liot-pliasc optical spectra compared with CMFGEN model with /3 = 4 (dotted line). Bottom panel shows the 
spectrum obtained with SCORPIO, and top and middle panels show the spectrum obtained with FOCAS. The model spectrum on the 
bottom panel is convolved with the 5A-wide Gaussian instrumental profile. 



photosphere of Romano's star coming to the two princi- 
pal conclusions. First, variability in this object is caused by 
correlated changes in mass-loss rate, wind velocity and hy- 
drostatic radius. Secondly, elementary abundances do not 
change significantly, we find similar helium and nitrogen 
overabundance in both states, characteristic for hydrogen- 
rich WNL stars, H/He~ 1.7-2.0 and N/He~ (3-5) x 10"^. 

We find that the bolometric luminosity of this object 
was higher during the eruption in 2005 by a factor of ■~1.5, 
that makes V532 one more example of an LBV that changes 
its luminosity. Together with the moderate intensity out- 
burst of AFGL2298, its behaviour indicates that even mod- 



erate amplitude LBV outbursts are accompanied by changes 
in bolometric luminosity. 
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Figure 6. Fitting the cool-phasc spectrum in the blue spectral range. Top panel: spectrum of V532 (Feb. 2005) (solid line) and our 
cool-phase model (convolved with the instrumental profile). Bottom: comparison of our cool-phase model spectrum with the archival 
spectra of the LBV star P Cyg. Spectra are normalized by the local continuum level. 



Table 4. Derived properties of V532 in the maximum and minimum brightness, and comparison with related stars in M33, Large 
Magellanic Cloud (LMC) and Milky Way (MW) galaxies, including the LBVs P Cyg and AG Car in visual minima (1990 December, 
2002 July). 



Star 


Gal. 


Sp. 


T, 






^2/3 


logL, 


log Afci 


f 




H/He 


Ref 






type 


[kK] 


[R©] 


[kK] 


[R©] 


[L©] 


[M0yr-i] 




[kin s-'-] 






WR124 


MW 


WN8h 


32.7 


18.0 






5.53 


-4.7 


0.1 


710 


0.7 


[1] 


WR40 


MW 


WN8h 


45.0 


10.6 






5.61 


-4.5 


0.1 


840 


0.75 


[2] 


WR16 


MW 


WN8h 


41.7 


12.3 






5.68 


-4.8 


0.1 


650 


1.2 


[2] 


R 84 


LMC 


WN9h 


28.5 


33.8 


24.9 


44.2 


5.83 


-4.40 




400 


2.5 


[4] 


BE 381 


LMC 


WN9h 


30.6 


20.8 


27.5 


26.0 


5.54 


-4.65 




375 


2. 


[4] 


AG Car 


MW 


WNll 


24.64 


67.4 


21.5 


88.5 


6.17 


-4.82 


0.1 


300 


2.3 


[5] 


Dec 1990 


























AG Car 


MW 


WNll 


18.7 


95.5 


16.4 


124.2 


6.0 


-4.33 


0.25 


195 




[5] 


Jul 2002 


























PCyg 


MW 


Blla+ 






18.7 


76.0 


5.8 


-4.63 


0.5 


185 


2.5 


[3] 


V532 


M33 


WN8 


34.0 


20.8 


31.7 


23.9 


5.7 


-4.72 


0.1 


360* 


1.9 




hot-phase 


























V532 


M33 


WNll 


22.0 


59.6 


20.4 


69.1 


5.89 


-4.4 


0.5 


200 


1.4 





cool-phase 



[1]- Crowther et al. (1999), [2]- Herald et al. (2001), [3]- Najarro (2001), [4]- Crowther & Smith (1997), [5] - Groh et al. (2009) 
* - foo was estimated using He I lines (Maryeva & Abolmasov 2010) 
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Figure 7. Fitting the cool-phasc spectrum in the red spectral range. Top panel: Spectrum of V532 (Feb. 2005) (solid line) and our 
convolved cool-phase model. Bottom: comparison of our cool-phase model with observation spectrum of LBV star P Cyg. Spectra are 
normalized by the local continuum level. 
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Figure 8. Approximation of the HeIlA4686 emission in the model 
spectrum by two Gaussians. The broad component of He IIA4686 
is shown by a thick line, narrow by a dashed line. 
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